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Intracellular transport is based on molecular motors that pull cargos 
along cytoskeletal filaments. One motor species always moves in 
one direction, e.g. conventional kinesin moves to the microtubule 
plus end, while cytoplasmic dynein moves to the microtubule minus 
end. However, many cellular cargos are observed to move bidi- 
rectionally, involving both plus-end and minus-end directed motors. 
The presumably simplest mechanism for such bidirectional transport 
is provided by a tug-of-war between the two motor species. This 
mechanism is studied theoretically using the load-dependent trans- 
port properties of individual motors as measured in single-molecule 
experiments. In contrast to previous expectations, such a tug-of- 
war is found to be highly cooperative and to exhibit seven different 
motility regimes depending on the precise values of the single motor 
parameters. The sensitivity of the transport process to small param- 
eter changes can be used by the cell to regulate its cargo traffic. 

Intracellular transport | motor regulation | cytoskeletal motors | stochastic 
processes 

The complex internal structure of biological cells depends 
to a large extent on targeted transport of vesicles, organelles 
and other types of cargo. This active intracellular transport 
displays the counterintuitive property that many cargos are 
observed to move bidirectionally, reversing direction every few 
seconds [1, 2]. This 'saltatory motion', which is faster and 
more persistent than Brownian motion, has been known for 
a long time [3]. With the improvement of experimental tech- 
niques, bidirectional motion has been found to be widespread, 
including particles such as mitochondria, pigment granules, 
endosomes, lipid-droplets, and viruses [2]. 

The long-range intracellular traffic inside biological cells 
is powered by molecular motors which transport cargos along 
microtubules (MTs). Some motors such as cytoplasmic dynein 
walk to the minus end, while others such as kinesin 1 or 2 walk 
to the plus end of the MTs. Cells often have a unidirectional 
MT cytoskeleton: The MT minus ends are typically located 
near the cell center, while the plus ends point outwards to- 
wards the cell periphery. Polarized cells like epithelial cells or 
axons possess a unipolar parallel MT array. Because of this 
unidirectional nature of the MT network and the motors, both 
plus and minus motors must be involved in the bidirectional 
transport of a single cargo. Indeed both kinesin and dynein 
are found simultaneously on various cellular cargos [4, 5, 6]. 
It is a matter of current research how the two motor species 
accomplish the bidirectional transport [1, 2, 7, 8, 9]. 

Two scenarios seem plausible [1, 2]: (i) Tug-of-war. Each 
motor species tries to move the cargo into its own direction, 
thereby performing a 'tug-of-war' on the cargo as depicted in 
Fig. 1. (ii) Coordination. An additional coordination com- 
plex prevents opposing motors from being active at the same 
time, thereby excluding state (0) in Fig. 1. In both cases, 
regulatory mechanisms, which may directly target the motors 
or the putative coordination complex, must be present in or- 
der to allow the cell to alter its motor transport in response 




Fig. 1. Cargo transport by 2 plus (blue) and 2 minus (yellow) motors: possible 
configurations (0), (+), and ( — ) of motors bound to the MT. For configuration (0), 
the motors block each other so that the cargo does not move. For configuration (+) 
and ( — ), the cargo exhibits fast plus and minus motion, respectively. 



to internal or external stimuli. 

The observed fast motion, and the complexity of bidirec- 
tional transport, as briefly reviewed in the following para- 
graphs, has led many authors to reject a tug-of-war scenario 
and search for a coordination complex. However, as shown in 
this article, this rejection of the tug-of-war scenario is prema- 
ture since a realistic tug-of-war leads to rather complex trans- 
port behavior that is not easily understood intuitively and, 
thus, may be erroneously interpreted as coordinated trans- 
port. 

Most quantitative data have been obtained experimentally 
in two model systems: pigment granule transport in fish and 
frog melanophores [10, 11, 9] and lipid-droplet transport in 
Drosophila embryos [12, 13, 14]. In melanophores, which are 
specialized pigment cells responsible for skin color, pigment 
granules move bidirectionally with similar velocities in both 
directions. They achieve net minus-end transport during an 
'aggregation period' because the average distance traveled in 
minus direction (the minus run length) is longer than the av- 
erage distance traveled in plus direction (the plus run length) . 
During a 'dispersion period', there is almost no net transport 
because of an increased minus run length [11]. 

Lipid-droplets are storage organelles for lipids. In late 
Drosophila embryos, they move on a unipolar MT array in 
the egg periphery. Their bidirectional motion exhibits differ- 
ent patterns in difi'erent stages of embryonic development. In 
particular, from so-called phase II to III their net transport 
direction changes from plus to minus because of an increase in 
the minus run length [12, 14]. This system is the only one for 
which force measurements have been performed so far. Stall 
forces have been found equal in plus and minus direction, in- 
dependent of the net direction of droplet transport [12, 14]. 

Various proteins that are necessary for the proper function 
or regulation of motor transport have been identified [15]. Ex- 
amples are the dynein cofactor dynactin, which is necessary 
for bidirectional transport in melanophores [16], or proteins 
like halo, klar or LDS2 in the lipid-droplet system [17]. 

Motor transport was found to be affected both by intra- 
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Fig. 2. A cargo with A^-|_ = 3 plus (blue) motors and A'^— = 2 minus (yellow) 
motors is pulled by a fluctuating number of motors bound to the MT. The configu- 
ration in the middle corresponds to (n-|-,n— ) = (2, 1). Only 5 out of 12 possible 
(n^, Tl— ) configurations are displayed. 



cellular regulation and by mutational changes in the motor 
structure. First, cellular regulation often leads to changes in 
only one direction. In the lipid-droplet system, net transport 
during embryogenesis is altered via a change in the plus run 
length [12], while in the melanophore system during skin color 
change the minus run length is changed [11]. Similarly, her- 
pesvirus capsids achieve targeting during entry and egress by 
modulation of the plus run length [18]. In all cases, the other 
direction is left unaltered. Second, mutation of the plus or mi- 
nus motor mostly causes reduced motion in both directions by 
decreasing run lengths or velocities, as observed by mutating 
dynein on lipid-droplets [13, 14] and kinesin on axonal pro- 
tein carrying vesicles [19]. However, in melanophores, kinesin 
inactivation leads to breakdown of plus motion and increased 
minus run lengths [11]. 

Interfering with the dynein-cofactor dynactin impairs 
transport in both directions in melanophores [16], but im- 
pairs minus and enhances plus transport of adenovirus parti- 
cles [20]. In the only m vitro experiment concerning bidirec- 
tional transport [21], a motility assay of kinesin and dynein, it 
was observed that increasing the number of dyneins enhances 
minus and impairs plus end transport. 

As shown in this article, all of these experimental observa- 
tions are consistent with the tug-of-war mechanism. In fact, 
we present the first explicit tug-of-war model that takes into 
account the experimentally known single motor properties and 
makes quantitative predictions for bidirectional transport. In 
our model, the motors act independently and interact only 
mechanically via their common cargo. We find seven possible 
motility regimes for cargo transport. Three of these regimes 
are dominated by the three configurations (0), {+), and ( — ) 
in Fig. 1 and represent no motion, fast plus motion, and fast 
minus motion of the cargo, respectively. The other motility 
states are combinations thereof; in particular, there are the 
two regimes ( — h) and (—0+) where the cargo displays fast 
bidirectional transport without and with pauses, respectively. 
During fast plus or minus motion, only one motor type is 
pulling most of the time and the tug-of-war appears to be 
coordinated. 

The different motility regimes are found for certain ranges 
of single-motor parameters such as stall force and MT affinity. 
Small changes in these parameters lead to drastic changes in 
cargo transport, e.g. from fast plus motion to bidirectional 
motion or no motion. We propose that cells could use the 



Table 1. Values of the single-motor parameters for ki- 
nesin 1, cytoplasmic dynein and an unknown plus motor kin? 
that transports Drosophila lipid-droplets. The kinesin 1 val- 
ues are taken from the cited references. The starred values 
are obtained by fitting experimental data of Drosophila lipid- 
droplet transport and agree with the cited references. 



parameter 


kinesin 1 


dynein 


kin? 


stall force Fs [pN] 


6 [29, 30] 


1.1* [12, 27] 
7 [31] 


1.1*[12] 


detachment force Fd [pN] 


3 [30] 


0.75* 


0.82* 


unbinding rate eo [s^^] 


1 [30, 32] 


0.27* [27, 33] 


0.26* 


binding rate ttq [s^^] 


5 [34] 


1.6* [33, 35] 


1.6* 


forward velocity vp [/im/s] 


1 [32, 36] 


0.65* [33, 37] 


0.55* 


back velocity vb [nm/s] 


6 [36] 


72* 


67* 



sensitivity of the transport to the single-motor properties to 
regulate its traffic in a very efficient manner. We illustrate 
this general proposal by providing an explicit and quantita- 
tive tug-of-war model for the lipid-droplet system. 

Results 

Model. To study the bidirectional transport of cargos, we de- 
veloped a model for a cargo to which A^+ plus and A''- minus 
motors are attached. Typically these numbers will be in the 
range of 1 to 10 motors as observed for many cargos in vivo 
[22, 12, 23]. For N+ = or N- = 0, we recover the model 
for cooperative transport by a single motor species as studied 
in [24]. We characterize each motor species by six parameters 
as measured in single molecule experiments (see Table 1 and 
text in the Supporting Information (SI)) as follows: it binds 
to a MT with the binding rate ttq and unbinds with the un- 
binding rate eo, which increases exponentially under external 
force, with the force scale given by the detachment force Fd. 
When bound to the MT, the motor walks forward with the 
velocity vp, which decreases with external force and reaches 
zero at the stall force Fs- Under superstall external forces, 
the motor walks backward slowly with backward velocity vb- 
The motors on the cargo bind to and unbind from a MT in 
a stochastic fashion, so that the cargo is pulled by n+ < N+ 
plus and n_ < A'^_ minus motors, where n+ and n_ fiuctuate 
with time, see Fig. 2. We have derived the rates for unbinding 
of one of the bound motors and for binding of an additional 
motor on the cargo from the single motor rates under the as- 
sumption that: (i) the presence of opposing motors induces 
a load force, and (ii) this load force is shared equally by the 
bound motors belonging to the same species (see SI text). 
We obtain a Master equation for the motor number proba- 
bility p{n+ , n_ ) that the cargo is pulled by n+ plus and n_ 
minus motors. The observable cargo motion is characterized 
by the motor states (n+ , n_ ) with high probability. If there 
is high probability for a state (714., 0) or (0, n_) with only one 
motor species bound, corresponding to Fig. 1(+) and ( — ), the 
cargo exhibits fast plus or minus motion, respectively. If there 
is high probability for a state with both motor species active, 
i.e. 71+ > and n_ > 0, the cargo displays only negligible 
motion into the direction of the motors that 'win' the tug- 
of-war, because the losing motors walk backward only very 
slowly. This corresponds to the blockade situation depicted 
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in Fig. 1(0). 

Motility states for the symmetric case. We first studied the 
instructive symmetric case, for which the number of plus and 
minus motors are the same and where plus and minus mo- 
tors have identical single-motor parameters except for their 
preferred direction of motion. Apart from being theoretically 
appealing, this symmetric situation can be realized in vitro if 
cargos are transported by a single motor species along antipar- 
allel MT bundles, and can also be used in vivo provided plus 
and minus end transport exhibit sufficiently similar transport 
characteristics. 

We solved our model for fixed motor numbers A'^+ = A*"- 
and fixed single-motor parameters and determined the prob- 
ability distribution p(n+,n_), see SI text. Depending on the 
values of these parameters, the model exhibits qualitatively 
different solutions, see Fig. 3, which we will call 'motility 
states' in the following. These motility states exhibit distinct 
cargo trajectories and velocity distributions as shown in Fig. 3 
and can formally be distinguished by the number of maxima 
of the motor number probability distribution p{n+,n-). This 
number of maxima is found to be either 1, 2, or 3. For the 
symmetric case, three types of maxima with the configura- 
tions of Fig. 1 occur: a maximum with only plus and no 
minus motors bound (+), one with only minus and no plus 
motors bound ( — ), and one with equal numbers of plus and 
minus motors (0). These maxima are found in the combi- 
nations (0), ( — h), and (— 0-f), leading to three qualitatively 
different motility states. 

(0) No motion For 'weak' motors with small stall to de- 
tachment force ratio / = Fs/Fd, the probability distribution 
p(n+, n_) has a single maximum at a state with an equal num- 
ber of bound plus and minus motors, see Fig. 3A1, and the 
velocity distribution has a peak at zero velocity, see Fig. 3A3. 
The corresponding cargo trajectories in Fig. 3A2 exhibit only 
small fluctuations around the initial position. This motility 
state (0) represents the blockade situation shown in Fig. 1(0) 
which one naively expects for a tug-of-war scenario. 

( — h) Fast plus and minus motion For strong mo- 
tors with large /, cargo movement is completely different. 
The cargo switches between fast plus-directed and minus- 
directed motion, see Fig. 3B2, and the probability distribu- 
tion p{n+,n-) has two maxima, see Fig. 3B1. At one maxi- 
mum only plus motors are bound to the MT (n+ > 0, n_ = 0) 
and at the other only minus motors (n+ = 0, n_ > 0), corre- 
sponding to the states (+) and ( — ) in Fig. 1 which are usually 
associated with coordinated transport rather than with a tug- 
of-war scenario. This behavior can be understood as follows: 
When more plus than minus motors are bound to the MT 
(n+ > n_), every plus motor experiences the force Fc/n+, 
while every minus motor experiences the larger force Fc/ri-, 
where Fc denotes the total force on the cargo. Since the un- 
binding rate increases strongly with increasing load force, mi- 
nus motors are more likely to unbind from the MT than plus 
motors, so that the predominance of the plus motors is fur- 
ther enhanced. After the unbinding of a minus motor the 
remaining minus motors experience an even larger force and 
are even more likely to unbind. As a consequence, the cargo 
experiences a cascade of minus motor unbinding events un- 
til no minus motor remains bound. A prerequisite for this 
unbinding cascade is that the motors can exert a sufficiently 
large force to pull off opposing motors from the MT, i.e. the 



stall force Fs has to be comparable or larger than the de- 
tachment force Fd- For small force ratios / = Fs/Fd, the 
pulling force has only a small effect on motor unbinding, so 
that no instability occurs and the cargo exhibits the blocked 
motility state (0). For large motor force ratio, the transient 
predominance of one motor type is thus amplified by a dy- 
namic instability and most of the time only one motor type 
is bound, as indicated in Fig. l(-f) and (— ). The emergence 
of cooperative behavior arising from the nonlinear force de- 
pendence of the unbinding rate has also been proposed as an 
explanation for collective effects in muscles [25] and mitotic 
spindle oscillations [26] . For the tug-of-war of 4 against 4 mo- 
tors with kinesin-like parameters as in Fig. 3B, about 90% 
of the time only one motor type is bound. During a plus or 
minus run, the effective velocity is however slightly reduced 
compared to the single-motor velocity by the sporadic binding 
and subsequent fast unbinding of an opposing minus motor. 
The velocity distribution in Fig. 3B3 has two peaks close to 
the single-motor velocities zbl/im/s. The direction of motion 
of the cargo is reversed when, due to a fluctuation, the de- 
feated motors become predominant. 

(—0+) Fast plus and minus motion with inter- 
spersed pauses Finally, in some intermediate parameter 
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Fig. 3. Motility states for the symmetric tug-of-war of N-\- = N— = 4 plus and 
minus motors. The three columns (A), (B) and (C) correspond to the three motility 
states (0), ( — h) and (—0+): (A) The no-motion motility state (0) is characterized 
by (Al) motor number probabilities p with a single maximum at an equal number 
of active plus and minus motors, {A2) trajectories with almost no motion, and (A3) 
velocity distributions with a single maximum at zero velocity (B) The motility state 
( — h) of fast bidirectional motion is characterized by (Bl) probabilities p with two 
maxima with only plus or only minus motors active, (B2) trajectories which exhibit 
switching between fast plus and minus motion, and (B3) bimodal velocity distribu- 
tions with two peaks close to the single-motor velocities of ib 1 pm/s. (C) The 
motility state (—0+) is characterized by (Cl) probabilities with three maxima corre- 
sponding to fast plus and minus and no motion, (C2) trajectories which exhibit fast 
plus and minus motion and pauses, and (C3) velocity distributions with three peaks. 
Both plus and minus motors in (B) have the kinesin 1 parameters of Table 1. The 
different motility behaviour in (A) and (C) is obtained by changing the single motor 
parameters in Table Ito (A) Fg = 2 pN, and (C) Fs = 4.75 pN, and €q = 0.4 s~^. 



The small peak near zero velocity corresponds to the no-motion states near the maximum for 
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are non-zero. 
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ranges, the probability distribution p{n+,n-) exhibits three 
maxima as shown in Fig. 3C1, a symmetric one corresponding 
to no motion as for motihty state (0) and two nonsymmetric 
ones corresponding to steady plus and minus motion as for 
state ( — +). The velocity distribution has three corresponding 
peaks, see Fig. 3C3, and cargo trajectories therefore exhibit 
bidirectional motion interspersed with pauses, see Fig. 3C2. 

Motility states for the asymmetric case. Bidirectional cargo 
transport in vivo is typically dependent on two different motor 
species for plus and minus motion. This plus-minus asymme- 
try can lead to net transport of the cargo in one direction. 
For example, in the motility state (—0-1-), the plus motion 
maximum (-I-) of the motor number probability can be larger 
than the minus motion maximum ( — ), see Fig. 4A1, which 
leads to longer plus runs compared to minus runs and to net 
plus motion of the cargo as illustrated by the trajectory in 
Fig. 4A2. The velocity distribution in Fig. 4A3 has the three 
peaks characteristic for the (— 0-f) regime but the peak at high 
plus velocity is larger than the one at high minus velocity. As 
cargo motion is no longer symmetric with respect to plus and 
minus motion, seven motility states are now possible, corre- 
sponding to the different combinations (+), ( — ), (0), ( — h), 
(0-I-), (-O)and (-0+) of the maxima {+), (-), and (0). The 
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Fig. 4. Motility states for the asymmetric tug-of-war of A'^-i- = 6 plus against 
A'^— = 6 different minus motors. The cargo is in one of seven motility states. 
The motility states (0), ( — h), and ( — 0-|-) are as for the symmetric case shown in 
Fig. 3, except that the plus-minus symmetry is lost as illustrated in column (A) for 
the ( — 0+) motility state. (A) The ( — 0-|-) motility state is characterized by (Al) 
three maxima in the motor number probability p at a plus, a minus and a no-motion 
state, (A2) trajectories with rapid plus and minus motion interspersed with pauses, 
and (A3) three peaks in the velocity distribution. Plus motion has a higher probability 
so that net motion is plus-end directed. (B) The (0-t-) motility state is characterized 
by (Bl) probabilities p with one maximum with only plus motors and one with plus 
and minus motors active, (B2) trajectories with fast plus motion and pauses, and 
(B3) velocity distributions with two peaks near the single plus motor velocity Vp^ 
= 0.55 flm/s and near zero. (C) The (-1-) motility state is characterized by (CI) 
probabilities p with a maximum with only plus motors active, (C2) trajectories with 
fast plus motion, and (C3) velocity distributions with a peak close to the single plus 
motor velocity. The motility states ( — 0) and ( — ) are similar to the states (O-j-) and 
(-1-) with plus and minus interchanged. (A) represents lipid-droplet transport: The 
plus and minus motors have the Drosophila plus motor (kin?) and dynein parame- 
ters of Table 1. The same parameters are used in (B) and (C) except for (B,C) Fs — 
= 0.45 pN and (B) eo- = 0-24 s^^ and (C) eo- = 54 s^^. 



new motility states (0-|-) and (-I-) are shown in Fig. 4B, C. 
The two other new states (—0) and (— ) are analogous with 
plus and minus motion interchanged. 

In the motility state (0-I-), the motor number probabil- 
ity has a maximum at the plus motion state (-f ) with only 
plus motors active and a maximum at the no-motion state (0) 
with both types of motors active, see Fig. 4B1. The corre- 
sponding velocity distribution in Fig. 4B3 has two peaks, one 
close to zero velocity and one at large plus motor velocity, and 
the cargo switches between fast plus motion and pauses, see 
Fig. 4B2. In the (-f) motility state in Fig. 4C, the motor 
number probability and velocity distribution exhibit a maxi- 
mum corresponding to fast plus motion^. 
In VIVO tug-of-war. To check whether our model can ac- 
count quantitatively for experimental observations, we applied 
our model to the bidirectional movements of lipid-droplets in 
Drosophila embryos. We chose this particular series of sophis- 
ticated experiments [12, 13, 14] because it is unique in pro- 
viding an estimate for the number of motors on the cargo, a 
high number of quantitative measurements of transport char- 
acteristics including cargo force measurements, as well as ob- 
servations in two different developmental phases (labeled wild 
type phase II and III, Wt II and Wt III) and in three different 
dynein mutation backgrounds. The droplets are transported 
by an unknown plus motor, presumably an unconventional 
kinesin, and cytoplasmic dynein [13]. 

We first considered the Wt II data. Cargo stall force mea- 
surements led to the conclusion that the droplets are on av- 
erage pulled by 5 plus and 5 minus motors, and that both 
types of motors have a single-motor stall force of 1.1 pN [12]. 
Since the number of active motors fluctuates stochastically, 
this should be the average number of pulling motors. We 
therefore fixed the total number of plus and minus motors to 
N+=N-= 6. 

We then performed simulations and varied the undeter- 
mined single motor parameters in order to fit the experimen- 
tally measured transport characteristics, namely plus and mi- 
nus run lengths, plus and minus stall forces, pause times after 
plus and minus travel, and plus and minus velocities of short 
and long runs, with an accuracy of ca. 10% (for the detailed 
procedure and results of this and the following fits see SI text 
and Table 2 and 3). The resulting parameters for dynein and 
the unknown plus motor (kin?) are listed in Table 1. The 
dynein parameters are in agreement with in vitro measure- 
ments of dynein properties when available. All other parame- 
ter lie in a reasonable range. The dynein backward velocity is 
an order of magnitude larger than for kinesin 1, in agreement 
with experiments [27, 28]. 

Fig. 4A shows a sample trajectory, the motor number 
probability and the velocity distribution for the droplet tug- 
of-war in Wt II. The cargo switches between fast plus and mi- 
nus motion and pauses but exhibits net plus motion because 
the probability for (-I-) states is higher than for (-) states. 
The cargo stall forces in plus and minus direction are equal, 
see SI Table 3. This shows that the cargo direction is not 
only determined by the motor forces but also by other motor 
properties, see Table 1. In this case, the higher plus motor 
detachment force makes it difficult to rip off the plus motors 
and thus favors plus motion. 



"This has also been found In the melanophore system [11]. 
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Fig. 5. Drosophila lipid-droplet transport in wild type phase M; tug-of-war of 6 Drosophila plus motors and 6 dyneins with parameters as in Table 1. (A) Distribution 
of run length, the distance traveled in one direction before a pause or a directional switch occurs. Minus (plus) run lengths are negative (positive). The gray bars are the run 
lengths observable with the experimental cutoffs of a minimum length of 0.16 s and 30 nm, while the white bars are obtained without the cutoffs and thus beyond experimental 
resolution. The lines are double exponential fits to the simulation data with decay lengths of ca. 0.1 flm and 1 ^m in both directions, of the same order of magnitude as in 
the experiments [14]. (B) Scatter plot of the run velocity (positive for plus and negative for minus runs) versus the run length of 500 runs in each direction shows a positive 
correlation: longer runs have higher velocities. The Spearman rank correlation coefficient both for plus and minus motion is larger than 0.7 with a significance level below 
10~ . Long runs have almost the maximal velocity which is the single motor velocity, Vjr^ in plus and — Vp— in minus direction. There are no data points for small run 
lengths and velocities because runs have been defined as periods with a velocity larger than 50 nm/s for at least 30 nm. (C) The correlation of run length and velocity can 
also be seen by considering short (0.03-0.1 /im) and long (0.5-1 /^m) plus end resp. minus end runs. Short runs have lower averages than long runs, which reproduces the 
experimental averages of [14]. 



A nontrivial consistency check of our model is provided by 
three additional features that we obtained from this model in 
close agreement with experimental observations even though 
these features were not used in order to determine the model 
parameters in Table 1. First, the distribution of plus and 
minus run lengths can be fitted by a sum of two exponen- 
tials^ with length scales of the same order of magnitude as 
obtained experimentally [14], see Fig. 5A. Second, the pause 
time distributions of pauses after plus and after minus runs 
are very similar and can be fitted by a single exponential func- 
tion with a time scale of the same order of magnitude as in the 
experiments [13], see SI Fig. 11. Third, there is a correlation 
between run length and run velocity: long runs have larger 
average velocities, see Fig. 5B. In the experiments [13, 14], 
this has been quantified by dividing the runs into short runs 
(0.03 — 0.1 /im) and long runs (0.5 — l^m). Short runs have 
approximately half the velocity of long runs, see Fig. 5C. In 
our model, this property refiects the correlation of the average 
number of active winning motors with the run length, see SI 
Fig. 8, and can be understood as follows. 

During a certain run, e.g. in the plus direction, minus mo- 
tors bind from time to time to the MT. This slows down the 
motion and causes a 'pause'. However, the active plus mo- 
tors generate a large force on this single minus motor, which 
is then ripped off fast from the MT. As a consequence, the 
'pauses' are too short to be detected experimentally and are 
only noticeable via the decreased average cargo velocity in the 
plus direction. If a cargo is pulled by many plus motors, this 
has two effects: (i) The effective cargo velocity is increased 
because opposing minus motors do not create large forces on 
each of the many plus motors and because the minus motor 
drops off very fast, (ii) The plus run length is larger because it 
is less probable that the minus motors take over. Both effects 
together lead to a correlation of run length and velocity. 
Mutation and regulation. Three different dynein mutations in 
Drosophila embryos of developmental phase II have been stud- 
ied in [13, 14], and all three lead to impairment of both plus 
and minus motion with reduced run lengths and stall forces. 
At first sight, this simultaneous impairment of both trans- 



port directions in response to mutations that affect only one 
of the two motor species seems to be inconsistent with a tug- 
of-war. However, using our tug-of-war model, we were again 
able to describe the observed behavior with an accuracy of 
about 10%. In order to do so, we only varied the minus motor 
paramaters and kept the plus motor parameters fixed to their 
Wt II values, see SI Table 3. 

In our model, the dynein mutations simultaneously modify 
several parameters of this motor, among which are its unbind- 
ing rate, its binding rate, and its detachment force. If only 
one of these parameters were modified, the resulting motor 
behavior would be easy to understand. First, if only the un- 
binding rate is increased, the minus motors unbind from the 
filament faster and thus produce less force on the plus motors, 
which leads to longer plus and shorter minus runs. Second, 
increasing only the minus motor binding rate has the oppo- 
site effect because dyneins are more likely to rebind to the 
filament. Third, if only the minus motor detachment force is 
enhanced, the ability of the minus motors to resist the plus 
motors is also enhanced, which increases minus and decreases 
plus run lengths. Therefore, if only a single parameter of the 
minus motor is modified, motion in one direction is enhanced 
whereas motion in the opposite direction is impaired. On the 
other hand, the overall effect of changes in several motor pa- 
rameters is difficult to anticipate intuitively and can lead to 
impairment of both directions as shown in our model. 

Furthermore, two different embryonic phases Wt II and 
Wt HI allow to assess the effect of cellular regulation. In 
Wt II net droplet transport is plus-end directed, while it is 
minus-end directed in Wt HI due to a reduction in plus run 
lengths. Apart from the stall forces, all other transport char- 
acteristics remain unchanged. We propose that the cellular 
regulation that causes this change targets the motor prop- 
erties. We therefore fitted the Wt HI data by varying the 
single-motor parameters as for the Wt II data. The fit shows 
that a tug-of-war can lead to impairment of motion in one 
direction while leaving the other direction unaffected, see SI 
Table 3. The obtained single motor parameters for Wt II and 
Wt HI are rather similar. This sensitivity of motion to the 



MiJIIer et al. 



single-motor parameters allows the cell to regulate its trans- 
port in an efficient way. 

In the only in vitro experiment concerning bidirectional 
transport [21], a motility assay of kinesin and dynein, it was 
observed that increasing the number of transporting dyneins 
enhances minus and impairs plus end transport. This is re- 
produced in our model, compare SI Table 3. 

Discussion 

We have investigated a model for bidirectional cargo trans- 
port based on a tug-of-war between plus and minus motors 
governed by mechanical interactions only. Our model exhibits 
many features that are usually attributed to a coordination 
mechanism. In particular, even for equally strong motors, a 
tug-of-war does not necessarily lead to the expected block- 
ade situation with almost no cargo motion as in Fig. 1(0), 
but can also lead to switching between fast plus end and fast 
minus end motion as in Fig. l(-l-) and (— ). This surprising 
behavior is caused by a dynamic instability arising from the 
strongly nonlinear force-dependence of the single-motor un- 
binding rate. This instability leads to a high probability of 
having only one motor type active at a given time. 

In our tug-of-war model, the motility behaviour of the 
cargo is very sensitive to the single motor properties. Chang- 
ing the motor stall force or MT affinity, for example, can 
lead to qualitatively different motility behaviour such as fast 
plus motion, no motion, or bidirectional transport. When 
we modified the single-motor properties to mimic the effect 
of either mutations or of regulatory processes, we found that 
motion in plus and minus direction can be affected in vari- 
ous ways. We found cases for which (i) motion was affected 



only in one direction, (ii) motion was impaired in one di- 
rection and enhanced in the other, and (iii) motion was en- 
hanced or impaired in both directions. This variability agrees 
with experimental observations where different systems also 
exhibit widely varying reactions to regulation or mutation 
[11, 12, 13, 14, 16, 18, 19, 20] as described in the introduction. 

Our tug-of-war model is thus in qualitative agreement with 
experimental data for bidirectional transport in vivo. Further- 
more, we have been able to quantitatively describe the exper- 
imental data for the Drosophila lipid-droplet system. The lat- 
ter system exhibits different transport regimes depending (i) 
on the different phases of the embryonic development, which 
are distinguished by distinct sets of regulatory proteins, and 
(ii) on the molecular structure of the motor proteins, which 
have been changed by mutations. In our theory, these dif- 
ferent transport regimes arise from variations in single motor 
parameters, but the basic transport mechanism underlying all 
of these regimes is provided by a tug-of-war between the two 
motor species. 

Our results show that the two scenarios for bidirectional 
transport displayed in Fig. 1, namely the tug-of-war and 
coordinated motor activity, are not mutually exclusive, but 
rather that the tug-of-war provides a mechanism for coordi- 
nated movement. 
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Model for a single motor 

A single motor can walk along a MT with velocity v, unbind 
from the MT with unbinding rate e and bind to the MT with 
binding rate tt. Our choice of the rates is based on the load- 
dependent transport parameters of single motors as measured 
in single-molecule experiments. When bound to the MT, the 
motor moves forward with a load-dependent velocity which 
decreases monotonically from the zero-load forward velocity 
vp to zero at the stall force Fa [1, 2, 3, 4]. For superstall 
loads the motor walks slowly backwards, as has been shown 
for kinesin [2, 5]. In our model, we used a piecewise linear 
force-velocity relation with 



v{F) 



VF (1 - F/Fs) 
vb{1-F/Fs) 



for 
for 



F <Fs 
F> Fs 



[1] 



Here, vb is the absolute value of the motor backward ve- 
locity. For forces smaller than the stall force, such a linear 
relation provides a good approximation for the experimen- 
tally determined force- velocity curves both for kinesin [1, 2, 3] 
and dynein [4]. For superstall forces, the shape of the force- 
velocity curve is not known precisely. In this range our linear 
relation can be considered as a Taylor expansion to first or- 
der in F — Fs- The detailed form of the force- velocity curve 
is however not crucial for our results, as long as it decreases 
monotonously and exhibits a small backward velocity. The 
unbinding rate of the motor from the MT increases exponen- 
tially with the applied force F: 



e{F)=eo exp(|F|/Fd 



[2] 



as measured for kinesin [1] and as follows from Kramers or Bell 
theory. The force scale is set by the detachment force Fd- The 
binding rate to the MT is difficult to assess experimentally. 
It is expected to depend only weakly on the load because an 
unbound motor relaxes and then binds from its relaxed state 
(see the discussion in [6]). We therefore take the binding rate 
equal to the zero- load binding rate ttq, independent of load: 



TV{F) = TTO 



[3] 



The single motor rates of kinesin 1, cytoplasmic dynein 
and an unknown plus motor that transports lipid-droplets 
in Drosophila embryos are shown in Table 1. For kinesin 1 
all parameters have been measured in single-molecule experi- 
ments, see the references in the table. For dynein, only part 
of the parameters have been measured, and for the stall force 
conflicting results have been reported by different labs, see 
the references given in Table 1. In addition, dynein is very 
sensitive to regulatory and accessory proteins [7]. The un- 
known dynein parameters and the parameters of the unknown 
Drosophila plus motor are obtained by fitting experimental 
data from Drosophila droplet transport [8, 9, 10] as described 
in the section 'Fit to the lipid-droplet data' below. 



Effective rates for the cargo 

We consider a cargo that is transported by constant numbers 
of N+ plus and A'"_ minus motors. At every time t, the state of 
the cargo with A'^+ plus and A'^_ minus motors firmly attached 
to it is fully characterized by the numbers 71+ and n_ of plus 
and minus motors that are bound to the MT and thus actively 
pull on the cargo at that time. The state of the cargo changes 
when a plus or a minus motor binds or unbinds to/from the 
MT, see Fig. 2. These changes are described by a Master 
equation for the probability distribution p{n+,n-,t) to have 
n-f- bound plus and n_ bound minus motors at time t. The 
rates of this Master equation describe the transitions corre- 
sponding to the arrows in Fig. 2 and are determined from the 
single-motor rates using the assumptions that the motors act 
independently and feel each other only due to two effects: (i) 
opposing motors act as load, and (ii) identical motors share 
this load. If each plus motor feels the load F+ (and generates 
the force —F+) and each minus motor feels the load —F- (and 
generates the force F-), this means that the force balance on 
a cargo pulled by n+ plus and n_ minus motors is 



n+F+ — ~n^F- 



Fc. 



[4] 



Here, the sign of the force is chosen positive if it is a load on 
the plus motors, i.e. if it points into the minus direction. If 
only one motor type is bound, i.e. if n+ = or n_ = 0, then 
F+ = F- = Fc = Q. A single bound plus motor thus feels the 
force F+ = Fc/n+. Using Eqs. [2] and [3], this implies that 
the effective rate for the unbinding of one plus motor is 

n+eo+ exp [FrJ{n+Fd+)] , [5] 

and the effective rate for the binding of one plus motor is 

(iV+-n+)7ro+. [6] 

Here and in the following, the index '-I-' labels plus motor 
properties. Analogous expressions hold for the minus motors 
with the parameters indexed by ' — '. 

The cargo force Fc is determined by the condition that 
plus motors, which experience the force Fc/tlj^, and minus 
motors, which experience the force —Fc/ri-, move with the 
same velocity, which is the cargo velocity Vc. 



Vc{n+,n-) =v+ {Fc/n+) = -V- {-Fc/n,-) 



[7] 



Here, the sign of the velocity is taken positive in the plus 
direction and negative in the minus direction. In order to 
have a unique solution Fc to this equation, both motors 
must have nonzero backward velocities; otherwise the single- 
motor force velocity relations v+{F) and V-{F) do not have 
well-defined inverses. In the case of 'stronger plus motors' 
n+Fa+ > n-Fa-, Eqs. [1] and [7] lead to the cargo force 
and velocity 



Fc{n+,n-) = Xn+Fa+ + {1- \)n-Fa. 

n+Fa+ ~ n-Fa- 
Vc(n+,n-) - 



n+Fa + /VF+ + n-Fa-/VB- 



with A = 1/ {1 + (n+Fa+VB-) / {n-Fa-VF+)). The cargo 
moves to the plus direction, Vc > 0. In the opposite case 
of 'stronger minus motors' with n+Fa+ < n-Fa~, in Eqs. [8] 
and [9] the plus motor forward velocity vf+ has to be re- 
placed by its backward velocity vb-, and the minus motor 
backward velocity vb- by its forward velocity vf- ■ The cargo 
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moves into the minus direction, Vc < 0. Typically the back- 
ward velocity is rather small, so that a cargo with n+, n_ > 
pulled by both types of motors usually moves very slowly, as 
indicated by the 'blockade' situation in Fig. 1(0). If however 
only one motor type is bound, e.g. if n_ = 0, the cargo moves 
fast with the single plus motor velocity Vc — Vf+ , correspond- 
ing to Fig. 1{+). 

Eqs. [5], [6], [8], and the corresponding equations for 
minus motors, fully determine the rates that enter the Mas- 
ter equation for the motor number probability p{n+,n- , t) on 
the state space < n+ < A''-f, < n_ < A''-. In each state 
{n+,n-), the cargo moves with velocity Uc(n+,n_) as given 
byEq. [9]. 

Experiments usually observe only cargos that have been 
bound to the MT for an unknown time period and monitor 
them over a timescale of minutes which is large compared 
to the times scales of motor (un-)binding, which are of the 
order of seconds. We are therefore interested in the long- 
time behaviour of the cargo which corresponds to the time- 
independent steady state solution p{n+,n-) of the Master 
equation. 

External forces 

If an external force fcxt is present, the force balance Eq. [4] 
becomes 



n+F+ 



-n-F- + Fc; 



[101 



Here again, forces are taken to be positive if they point into 
the minus direction. Carrying through the same calculation 
as for the case without external force leads to the velocity of 
a cargo transported into the plus direction by n+ active plus 
and n_ active minus motors under an opposing external force 
-Foxt: 



Wc(n+,n_,Fc: 



n+Fs. 



■n-Fs. 



fc, 



n+Fs+/vF+ + n-Fs-/vB- 



[11] 



In the case of minus motion under an opposing force Foxt 
(which is then negative), the plus motor forward velocity vf+ 
has to be replaced by its backward velocity vb-, and the 
minus motor backward velocity vb- by its forward velocity 
vf--^ Thus the stall force of a cargo pulled by n+ plus and 
n_ minus motors is given by 

Fsc{n+,n-) = n+Fs+-n-Fs-, [12] 

as intuitively expected. 

Numerical calculations 

For a given set of single-motor parameters, we numerically 
calculate the steady state solution p{n+,n-) as the nuUspace 
of the transition matrix of the Master equation [11]. Such 
motor number probabilities are shown in Fig. 3A2-C2 and 
Fig. 4A2-C2. We then determine the locations (n+,n_) of 
maxima of p(n+, n_) which define the 'motility state' of the 
cargo. p{n+ , h- ) is a maximum if it is larger than its direct 
and diagonal neighbours on the state space {n^,n-). A maxi- 
mum at a state (n+ , 0) with n_|. > corresponds to a plus mo- 
tion state labeled by '-|-', a maximum at (0, n_) with n_ > 
to a minus motion state labeled by ' — ', and a maximum at 
a state {h+,h-) with both n+ and n_ larger than zero to a 
no- motion state labeled by '0'. For a given set of single-motor 



parameters, there is at most one of each type of maximum. 
Thus there are seven possible combinations of maxima, which 
give the seven motility states (-I-), ( — ), (0), (0+), (—0), ( — h), 
and (—0-1-). In the motility state (0-I-), for example, there are 
two maxima, one at a plus-motion state (n+,0), and one at a 
no-motion state {n+,n-). If the probability maximum is at 
(0, 0), the cargo is in the 'unbound' state. 

In order to obtain dynamical quantities such as trajecto- 
ries, run lengths and run velocities, we generate individual 
cargo trajectories using the Gillespie algorithm [12] for the 
binding/ unbinding dynamics and let the cargo move with the 
velocity Vc in the intervals between binding/unbinding events. 
At the start of the simulation, a cargo is bound to the fila- 
ment by a random number of n+ active plus and n_ active 
minus motors. In order to suppress transient behavior due to 
initialization bias, measurement of run lengths and velocities 
is started only after a time lapse of lO* s. The obtained ve- 
locity and run length distributions as shown in Fig. 3A3-C3, 
Fig. 4A3-C3 and Fig. 5A, and all average values, are obtained 
from 20-50 trajectories, each of which lasts 10'' s. 

For comparison with experiments, the analysis of the sim- 
ulated trajectories is performed in close analogy with the anal- 
ysis of the experimental trajectories of [8, 9, 10]. The cargo 
displacement (as shown in Fig. 3A2-C2 and Fig. 4A2-C2) is 
recorded at video frequency of 30/s. The obtained trajectories 
are segmented into runs and pauses by using the definitions 
from [9, 10]: A cargo is considered to be moving into the 
plus (minus) direction if its velocity Vc is larger than 50 nm/s 
(smaller than —50 nm/s) and pausing else. A run has to be 
at least 30 nm and 0.16 s long, and a pause must be longer 
than 0.23 s and cover a distance smaller than 30 nm. The run 
velocity is defined as the ratio of run length and run time. 

Motility diagrams 

In the symmetric tug-of-war, for which the number of plus 
and minus motors are the same and for which plus and mi- 
nus motors have identical single-motor parameters except for 
their preferred direction of motion, the cargo motion depends 
on four dimensionless parameters: the number N = 7V+ — N- 
of plus and minus motors, the stall force to detachment force 
ratio / = Fs/Fd, the MT desorption constant K = eo/vro, and 
the backward- forward velocity ratio v — vb/vf- Depending 
on the values of these parameters, the cargo is in one of three 
possible motility states: (0), ( — h), and (—0-1-). These states 
are characterized by qualitatively different motility behaviors 
as illustrated in Fig. 3 for three sets of parameters. 

Regulation of cargo motion requires change of the motor 
parameters. To show the effect of parameter changes explic- 
itly, we calculate the 'motility diagram' shown in Fig. 6 for 
the symmetric tug-of-war with TV = N+ — N- = 4 sym- 
metric plus and minus motors. This is done as follows: The 
single-motor parameters are taken to be equal to the kinesin 1 
values as given in Table 1 except for the plus and minus 
motor unbinding rates eo = ^o+ = eo- and the stall forces 
Fs = Fs+ = Fs~. The parameter space {eo,Fs) is then 
explored systematically'^, and for each point the maxima of 



The case of assisting force is not treated here. It needs a definition of the single motor force- velocity 
relation Eq. [ 1 ] also for assisting loads i^ < 0. 

All other single-motor parameters are kept constant. However, as far as backward motion is con- 
cerned, not the backward velocity v^ is kept constant but rather the backward slope v^/Fs of 
the force-velocity curve, with F's equal to the kinesin 1 value. 
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the motor number probability p{n+,n-) is calculated and the 
cargo motility state is determined, as described in the section 
'Numerical calculations' above. When the maxima between 
two scanned points change, we zoom in between these points 
in order to determine the transition point more accurately. 
The lines shown in Fig. 6 consist of these points at which the 
locations of the maxima change. 

For large desorption constants K, the motors have a low 
affinity to the MT; therefore the number of bound motors 
in Fig. 6 is low for low K. For very high desorption con- 
stants K larger than the number of motors A'' = 4, the cargo 
is 'unbound', i.e. the maximum of the motor number prob- 
ability is at {h+,h-) = (0,0). For small force ratios /, the 
probability distribution p(n+,n_) has a single maximum at 
a state (n, h) with an equal number n = n+ = n_ of active 
plus and minus motors and is in the no-motion motility state 
(0) (green). For large force ratios /, the motors can generate 
forces large enough to rip off opposing motors since the stall 
force is large compared to the detachment force. This leads 
to the unbinding cascade described in the main text, and the 
motor number probability has maxima at states with only one 
active motor type, i.e. at (n+, 0) with only active plus motors 
and at (0, n„ ) with only active minus motors. In the latter 
situation, the cargo is in the ( — h) motility state (yellow). For 
intermediate values of /, both types of maxima coexist, and 
the cargo is in the (—0+) motility state (red). 

The lines in Fig. 6 that separate regions of different color 
mark the parameter sets where maxima of the motor num- 
ber appear or disappear and correspond to bifurcation lines. 
For example, when passing from the cross labeled A in the 
(0) regime with a maximum at (n+, n_) = (3, 3), correspond- 
ing to no motion, into the (—0+) regime by increasing the 
force ratio /, two new maxima at {h+,n-) = (4,0) and 
{n+,n-) = (0,4) appear. This means that the cargo now ex- 
hibits fast plus and minus motion. However, this fast motion is 
rarely observed since the two new maxima have less probabil- 
ity than the pause-maximum at (3,3). With further increase 
of /, the maxima at (4,0) and (0,4) grow, while the maximum 
at (3,3) shifts to (2,2) and diminishes. Cargo motion becomes 
more and more dominated by fast plus and minus motion, 
and pauses become less frequent until they finally disappear 
together with the maximum at (2,2) when passing the line to 
the ( — +) regime. At the cross labeled C in Fig. 6, the cargo 
switches between fast plus and minus motion. A direct tran- 
sition from the (0) to the { —+) regime (without passing the 
(—0-1-) regime) occurs only for maxima with low motor num- 
bers, when the no-motion maximum and the two fast motion 
maxima are neighbours (either direct or diagonal) on the dis- 
crete state space {n+,n^). This happens, for example, when 
passing from a region with a maximum at {n,n) = (1, 1) to 
a region with maxima at (n, 0) = (1,0) and (0, n) = (0,1). 
The bifurcation line then marks the parameter set where all 
maxima have the same probability. 

We investigated the average times between direction rever- 
sals in the ( — +) motility state and found it to grow exponen- 
tially with the motor number. This indicates that the tran- 
sitions between the different motility states become nonequi- 
librium phase transitions in the limit of large motor numbers. 

Fig. 7 shows the motility diagram for the tug-of-war of 
7V-I- = 6 plus against A''_ = 6 minus motors with parame- 
ters closely related to the Drosophila lipid-droplet transport 



in wild type phase II. The diagram does not show maxima 
locations but only the bifurcation lines between the different 
motility states. It is generated in a similar way as Fig. 6: 
The single-motor parameters correspond to the values of the 
Drosophila plus motor kin? and dynein as given in Table 1 
except for the minus motor unbinding rate eo- and stall force 
Fs-. The parameter space (eo-,-Fs_) is then explored sys- 
tematically in the same way as for the symmetric case. Fig. 7 
exhibits seven motility states. The 'new' motility states (-f), 
( — ), (0+) and (—0) are asymmetric with respect to plus and 
minus motors and, thus, were not present for the symmetric 
tug-of-war. Fig. 7 shows that minus motors with high afHn- 
ity to the MT, i.e. with low desorption constant K-, favor 
minus motion. A high force ratio /_ enhances the unbinding 
cascade that leads to fast motion in the plus and/or minus 
direction. Minus motors with high MT affinity but low force 
ratio tend to block motion and lead to pauses. The irregular 
shape of the bifurcation lines between the motility states is 
a discretization effect and corresponds to transitions between 
different locations of the maxima of the motor number proba- 
bility p, similar to the changes of maxima locations shown in 
Fig. 6. 



Fit to the lipid-droplet data 

We applied our tug-of-war model to the bidirectional trans- 
port of lipid-droplets in Drosophila embryos as studied exper- 
imentally in [8, 9, 10]. Various transport characteristics were 
measured in two different embryonic phases (labeled wild type 
phase II and III, Wt II and Wt III) and three different dynein 
mutation backgrounds (labeled -D/ic®~^°/+, Dhc^~^ /+ and 
Dhc^-^/Dhc^-^°) during phase II. 

We first considered the Wt II data. From cargo stall force 
measurements, the experimenters concluded that the droplets 
are pulled in the plus and minus direction by five plus and 
five minus motors, respectively. Since the number of active 
motor fluctuates stochastically, these numbers represent the 
average number of pulling motors. We therefore chose the to- 
tal number of plus and minus motors to be A'^+ — N- — 6. 
The droplets are transported by dynein [9] and an unknown 
plus motor which we call kin?. The single motor parameters 
of dynein are only partly known, and for the stall force dif- 
ferent labs have reported different results, see Table 1. In the 
droplet experiments, cargo stall force measurements indicate 
a single motor stall force of 1.1 pN for both plus and minus 
motors"^ [8]. We used this value for the stall forces of both 
motors and varied the remaining 10 single motor parameters 
Fd±, eo±, 7i"o±, fF± and vb± in order to fit the experimental 
data. 

We generated and analyzed cargo trajectories as described 
above in the section 'Numerical calculations'. In particular, 
we used the experimental time resolution* and definition of 
runs and pauses, including the experimental cutoffs. The 
choice of the time resolution and the cutoffs is crucial since 
they strongly affect the magnitude of average run lengths, ve- 
locities or pauses. For example, short runs or pauses, which 
are easily accessible in simulations, may be unobservable in 



■ For dynein, this value Is In agreement with the stall force reported by [13]. The low stall force 
for the unknown plus motor Implies that this motor should be different from kinesin 1 because the 
kinesin 1 stall force is 6 pN. 
The spatial resolution in the experiments Is of the order of nanometers and therefore unproblematic. 
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experiments because they are below time resolution or buried 
in noise. For the experimental determination of plus and mi- 
nus droplet stall forces, the droplets had to be moving in a 
given direction 'for a few seconds' in order to decide their 
direction before the measurement was performed [8, 9]. We 
therefore averaged stall force values, calculated according to 
Eq. [12], only over 'very long' runs that last more than 3 
s. The experiments distinguished pauses after plus and after 
minus runs. We adopted this distinction although both in 
the experiments and in our simulations both types of pauses 
are very similar, see below. Furthermore, the experimenters 
defined 'short runs' of length 30-100 nm and 'long runs' of 
length 500-1000 nm, and calculated average velocities of both 
types of runs. We followed this procedure. 

For fitting, we compared 10 transport characteristics, 
namely plus and minus run lengths, plus and minus stall 
forces, pause times after plus and minus travel, and plus and 
minus velocities of short and long runs^ as shown in Table 3. 
For this purpose, we defined a 'distance function' between 
model and experiment as the sum of squared differences be- 
tween the experimentally measured and simulated quantities. 
As the different quantities are of different order of magnitude, 
they were rescaled in such a way that the experimental values 
are of order unity. We then minimized this distance function 
with respect to the unknown model parameters. 

For the Wt II fit, these are the 10 unknown single mo- 
tor parameters listed above. We first chose 'reasonable' set 
of parameters. Here 'reasonable' means that the motor pa- 
rameters must be of the order of magnitude of experimental 
single-motor parameters and that the simulation results must 
be of the order of magnitude of the experimental results. We 
then used the Nelder-Mead downhill simplex algorithm [14] 
to minimize the distance function starting from this initial 
choice. As this is only a local minimum, we repeated this 
procedure for several starting parameter sets until we found a 
minimum that reproduces the experimental data within about 
10%. 

In wild type phase III (Wt III) , reduced stall forces led to 
the conclusion that the average number of motors pulling the 
cargo in both directions is only four motors. We therefore took 
the total number of motors on the cargo to be A*'-|- = A''- = 5. 
For the Wt III fit, the single-motor parameters were set to the 
Wt II values. Then the simplex algorithm was started from 
this parameter set to minimize the 'distance function' with 
the Wt III experimental values. 

For the dynein mutation fits, only the six minus motor 
parameters Fs-, Fd-, eo-, 7i"o-, vf-, vb- were used as fitting 
parameters because the mutation only affects the dyneins, and 
does so in an unknown way. The plus motor parameters were 
kept fixed at the values from the Wt II fit since droplet mo- 
tion in the dynein mutation background were investigated in 
embryonic phase II. 

The single motor parameters resulting from all these fits 
are shown in Table 2 (the Wt II values are also shown in 
Table 1). They are within the expected range of motor pa- 
rameters. The unbinding and binding rates are of the order 
of 1 s~^ as measured for motors like kinesin 1 [15, 16, 17] 
and kinesin 3 [18]. For dynein, the unbinding and binding 
rates lie in the experimental range [19, 13, 20]. The forward 
velocities are of the order of 0.5 ^m/s which is close to the 
droplet velocity measured during long runs. This means that 



the tug-of-war does not substantially reduce the single-motor 
velocity. The backward velocity is two orders of magnitude 
smaller than the forward velocity but one order of magnitude 
larger than the kinesin 1 backward velocity. For dynein, this 
is in agreement with experiments [13, 21]. The wild type de- 
tachment forces obtained from the fit are approximately half 
of the stall force, similar as for kinesin 1, for which however 
both force scales are larger, see Table 1. 

A comparison of the experimental data and the corre- 
sponding fit result simulation data is shown in Table 3. They 
all agree within 10%. 

A remarkable feature of droplet transport is the positive 
correlation of run length and velocity, see Fig. 5B: longer runs 
have larger velocities. The correlation persists when consider- 
ing run times and velocities instead, see Fig. 8A. This is more 
meaningful since run length and velocity are trivially linearly 
correlated due to the fact that high velocities lead to larger 
displacement. As explained in the main text, the correlation 
is caused via a correlation of the average number of active 
winning motors with the run length, see Fig. 8B. In the ex- 
periments, the correlation has been quantified by dividing the 
runs into short and long runs as defined above and comparing 
the average velocities of these runs. For visual comparison, 
both model and experimental results are shown in Fig. 9. 

Although not used in the fitting procedure, the distribu- 
tions of run length and of pause time show the same qual- 
itative and similar quantitative behavior in simulation and 
experiment. Fig. 5A and Fig. 10 show the plus and minus 
run length distributions for all genotypes in wild type phase 
II, which can all be fitted by a double exponential function. 
The same behavior was found in the experiments [10]. The 
short and long decay lengths of these fits are listed in Table 4. 
The short decay lengths are ca. 0.1 fim, while the long decay 
lengths are of the order of 1 fj,in and vary in the different ge- 
netic backgrounds. Although not used for fitting, simulation 
and experimental values are of the same order of magnitude 
and agree within 50%. 

Fig. 11 shows the pause time distribution for Wt II pa- 
rameters. We did not distinguish pauses after plus and minus 
runs here because they were statistically indistinguishable. 
The pause time distribution can be fitted by a single expo- 
nential distribution with the time scale 0.38 s. All this is also 
found for the experimental distributions of pauses after plus 
and after minus runs, which are very similar to each other and 
can be fitted with a single exponential function with the time 
scales 0.24 s and 0.29 s, respectively [9]. 

The agreement of simulation and experiment shows that 
our tug-of-war model can describe the lipid-droplet data. Two 
features are particularly remarkable because they are not ex- 
pected within a naive picture of a tug-of-war: First, the Wt II 
and Wt III data represent motion with balanced stall forces 
in plus and minus direction, but the motion is net plus end- 
directed in phase II and net minus-end directed in phase III. 
Second, the dynein mutation data, which exhibit an impair- 
ment of both plus and minus motion, could be reproduced by 
varying the dynein single-motor parameters only. Both obser- 



"^We did not use quantities that are extremely sensitive to the detectability of pauses, such as 
the percentage of direction reversals associated vyith pauses, or the average time betvyeen pauses, 
because pauses in simulation are more readily detected than in experiment. We also did not use 
'secondary' quantities that were obtained by further processing of the data, such as fits to run length 
or pause time distributions or quantities calculated from these fits. 
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vations are in agreement with our tug-of-war model, as shown 
by the successful fit. The reason is that in our model cargo 
motion depends on six different motor properties (stall and 
detachment force, binding and unbinding rate, and forward 
and backward velocity) for each motor type, which leads to a 
rather complex behavior. In particular, cargo motion is not 
only determined by the motor forces but also by other motor 
properties, which leads to a variable response to perturba- 



tions such as mutation or regulation. As shown in Table 3, it 
is possible (i) to change only one direction and leave the other 
direction unaffected (Wt III, or change of only minus motor 
stall force or forward or backward velocity) (ii) to impair one 
direction and enhance the other (change of only minus mo- 
tor unbinding or binding rate or detachment force), or (iii) to 
impair both directions (dynein mutations). 
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Fig. 6. Motility diagram for the tug-of-war of A^_j_ = N— = 4 symmetric plus and minus motors with identical single-motor parameters but different forward directions. 
Depending on the motor force ratio / = Fs/^d of the stall and the detachment force and the MT desorption constant K = eo/TTo, the cargo transport exhibits three 
different motility regimes denoted by (0), ( — h), and (—0+). These regimes are defined via the number and locations of the maxima of the motor number probability 
distribution p(n_|_, ri_) as described in the text. The lines in the motility diagram separate regions in which the maxima of the probability distribution are located at different 
motor number states. The colors separate regions with different motility states. In motility state (0) (green) the motor number distribution has a single maximum at a 
no-motion state with an equal number of plus and minus motors bound at {h,fl) with 1 < ri < 4. The two neardiagonal maxima at (1,2) and (2,1) are also counted as a 
single diagonal maximum, which in a continuous state space would be at (n, n) with 1 < 72 < 2. If the maximum is at (n, n)=(0,0) the cargo is considered as 'unbound' 
(gray). In the ( — \-) regime (yellow), the probability distribution exhibits two maxima with only plus or only minus motors bound at (72,0) and (0,77.) with 1 < n < 4. 
The cargo is in the (—0+) regime (red) if the probability distribution exhibits three maxima. The parameters are kinesin-like as in Table 1 except for the stall force Fs and 
the unbinding rate eo, which are varied. The crosses labeled A, B and C mark the parameter sets for the cargo trajectories, probability and velocity distributions shown in the 
corresponding Fig. 3A-C, which illustrate the qualitatively different motility behaviours for the different motility regimes. The cross B in the (- — h) regime corresponds to the 
complete set of kinesin parameters with / = 6/3 and K = 1/5, while A in the (0) regime is at / ^ 2/3 and K = 1/5 and C in the ( — 0-f) regime is at 4.75/3 and 
K = 0.4/5. 
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Fig. 7. Motility diagram for the asymmetric tug-of-war of 7V+ ^ 6 plus against 7V_ ^ 6 minus motors related to the Drosophila lipid-droplet transport. Plus and 
minus motor parameters are for the Drosophila plus motor (kin?) and dynein, respectively, as in Table 1 except for the stall force Fs— ^nd the detachment rate eQ_ of the 
minus motors which are varied. Depending on the minus motor force ratio /_ = Fs — / F^— and the MT desorption constant K— ^= eo— /ttq— , the cargo exhibits seven 
different motility regimes: no motion (0), fast plus motion ( + ), fast minus motion ( — ), coexistence between two of these states as indicated by (0+), (—0), and ( — h), 
and three-state coexistence (—0+). The motility states are defined via the maxima of the motor number probability distribution p{Ti-\- , n_). Plus motion (+) corresponds 
to a maximum at (n-^,0) with only plus motors active, minus motion ( — ) to a maximum at (0,72—) with only minus motors active, and no motion (0) to a maximum at 
(n_^,fi_) with both motor types active. The other motility states exhibit the combinations of these maxima as indicated by the notation. For example, in the regime (O-h), 
the probability distribution has one maximum at a plus motion state and one maximum at a no-motion state. The crosses labeled A, B and C correspond to the parameter values 
for Fig. 4A-C, which illustrate the qualitatively distinct motility behavior for the different motility regimes. The cross A in the (—0+) regime corresponds to the tug-of-war for 
Drosophila droplet transport in wild type phase II with parameters as given in Table 1. The cross B in the (0+) regime is at /— ^ 0.60 and K— = 0.15, and C in the 
(+) regime at /_ = 0.60 and K- = 0.34. 
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Fig. 8. Scatter plots (A) of run velocities against run times and (B) of average number of active plus resp. minus motors against absolute run lengths for the Drosophila 
droplet tug-of-w/ar in wild type phase II for 500 runs in each direction, divided up into the positive plus runs (red) and the negative minus runs (blue). (A) The absolute run 
velocities are larger for longer runs and almost reach their maximal values of the single motor velocities Vp^ = 0.55 /im/s and Vp^ = 0.65 /Jm/s for very long runs. This 
show/s that the correlation of run length and velocity discussed in the main text persists when considering run times and velocities instead. This is more meaningful since run 
length and velocity are trivially linearly correlated due to the fact that high velocities lead to larger displacement. (B) The reason for the correlation is that longer runs also 
have a higher average number of active pulling motors, compare the discussion in the main text. There are no data points for small run times, lengths and velocities because 
runs have been defined as having a velocity of at least 50 nm/s for at least 30 nm and 0.16 s. 
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Fig. 9. Mean run velocity for short (30-100 nm) and long (500-1000 nm) runs for different Drosophila embryonic phases (Wt II and Wt III) and dynein mutations 
(Dhc^^ /+, Dhc /+ and Dhc / Dhc ) for (A) plus-end motion and (B) minus-end motion. Short runs have shorter velocities than long runs, and changes 

for the different phases or genotypes are via the velocity of the long runs. The simulations values (black) are taken from Table 3 and agree well with the experimental values 
(gray), which are read off from Fig. 4 in Ref. [14] and Fig. 7 in Ref. [13]. 
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Fig. 10. Run length distributions for three different dynein mutations labeled (A) Dhc^''^^ /+. (B) Dhc^^'^ /+, and (C) Dhc^^''^ / Dhc'^^^° . Plus (minus) run 
lengths are positive (negative). In all cases, the run length distribution can be fitted by a double exponential function corresponding to the solid line. The single motor 
parameters of the corresponding simulations are given in Table 2. The two decay lengths of the double-exponential fits are listed in Table 4. 
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Fig. 11. The pause time distribution (gray Inistogram) for tine tug-of-war with wild type phase II parameters can be fitted by a single exponential distribution (line) with 
time scale 0.38 s. We did not distinguish pauses after plus and after minus motion because they are statistically indistinguishable. The full graph is a semi-logarithmic plot 
of the original distribution shown in the inset. The first bar has reduced counts because of the experimental definition of a pause to be longer than 0.23 s. The distribution 
looks similar to the experimental distributions of pauses after plus and after minus motion in Fig. 5 of Ref. [13], which have been fitted with single-exponential distributions 
with time scales 0.24 s and 0.29 s. These time time scales are smaller than the experimental average pause times of 0.55 s for pauses after plus and 0.62 s for pauses after 
minus motion. Similarly, in simulation, the time scale 0.39 s of the exponential fit is smaller than the average pause time of 0.61 s. This indicates that the distribution is in 
fact not single exponential but has another, smaller time scale. This smaller time scale is below experimental resolution and shows up only when the experimental lower time 
cutoff 0.23 s for the pauses is removed (white histogram). 



Table 2. Single-motor parameters for the fits to Drosophila lipid-droplet transport in wild type phase II (Wt II) and III 
(Wt III), and for three different dynein mutations {Dhc^^^°/+, L>ftc* V+' ^"d Dhc^'^ / Dhc'^'^°). The motor numbers are 
A*'-|^ — N- — 6 except for Wt III with N+ = A'^ = 5. The plus motor parameters for the dynein mutations are as for the plus 
motors in Wt II. 





Wt II 




Dhc''-^''/+ 


Dhc'-'/+ 


Dhc'-^/Dhc''-^" 


Wt III 




Motor direction 


plus 


minus 


minus 


minus 


minus 


plus 


minus 


stall force Fs [pN] 


1.1 


1.1 


1.1 


1.1 


0.85 


1.1 


1.1 


detachment force Fd [pN] 


0.82 


0.75 


0.88 


0.84 


1.1 


0.82 


0.81 


unbinding rate eo [s~^] 


0.26 


0.27 


0.45 


0.37 


0.54 


0.25 


0.27 


binding rate ttq [s^^] 


1.5 


1.6 


1.8 


2.0 


1.8 


1.4 


1.6 


forward velocity vp [^im/s] 


0.55 


0.65 


0.69 


0.49 


0.44 


0.55 


0.63 


back velocity vb [nm/s] 


67 


72 


77 


76 


53 


68 


73 
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Table 3. Mutation and regulation in lipid-droplet transport: results of the fit to the Drosophila lipid-droplet data. The 
first 10 lines show a comparison of the average plus and minus stall forces, the average plus and minus run lengths, the 
average times of pauses after plus and after minus runs, and the average velocities of short and long runs in plus and in 
minus direction as obtained in simulation (sim.) and experiment* (exp.) for wild type phase II (Wt II) and III (Wt III), 
and for three different dynein mutations {Dhc^^^° /+, Dhc^^^ /+, and Dhc^^^ / Dhc^^^"). The last column describes the net 
effect on motion (run lengths and velocities) as compared to the Wt II values. The last seven lines show the effect of a 
change of one single motor parameter from the Wt II value, given in front of the arrow, to the value given after the arrow 
in the first column. 
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effect 






pN] 


[/^ 


m] 
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s] 




[m 


Tl/S] 




on 












after 


short 


runs 


long 


runs 


motion 




plus 


minus 


plus 


minus 


plus 


minus 


plus 


minus 


plus 


minus 


plus/minus 


Wt II (sim.) 


5.4 


5.3 


0.84 


0.66 


0.51 


0.61 


0.21 


0.21 


0.42 


0.47 


-/- 


Wt II (exp.) 


5.5 


5.5 


0.84 


0.62 


0.55 


0.62 


0.21 


0.22 


0.43 


0.47 




D/ic''-i7+ (sim.) 


5.2 


5.0 


0.58 


0.53 


0.54 


0.54 


0.21 


0.23 


0.41 


0.50 


impaired/impaired 


D/ic«-"7+ (exp.) 


5.5 


5.5 


0.56 


0.49 


0.60 


0.62 


0.21 


0.22 


0.46 


0.55 




D/ic*^"V+ (sim.) 


5.3 


5.1 


0.41 


0.32 


0.66 


0.66 


0.20 


0.21 


0.39 


0.37 


impaired/impaired 


D/ic*-V+ (exp.) 


5.1 


5.5 


0.38 


0.29 


0.71 


0.70 


0.21 


0.23 


0.37 


0.40 




Dhc^-^/Dhc''-^^ (sim.) 


5.0 


3.9 


0.29 


0.15 


0.71 


0.75 


0.14 


0.21 


0.29 


0.33 


impaired/impaired 


Dhc^-^/Dhc''-^° (exp.) 


4.7 


3.7 


0.31 


0.17 


0.71 


0.76 


0.17 


0.20 


0.29 


0.33 




Wt III (sim.) 


4.3 


4.4 


0.42 


0.60 


0.59 


0.59 


0.20 


0.20 


0.40 


0.45 


impaired/- 


Wt III (exp.) 


4.4 


4.4 


0.42 


0.60 


- 


0.60 


- 


0.20 


- 


0.44 




Fs- = l.lpN ^0.8pN 


5.3 


3.9 


0.83 


0.25 


0.75 


0.83 


0.16 


0.24 


0.39 


0.47 


-/impaired 


Fd-= 0.75 pH ^l.OpN 


4.9 


5.3 


0.24 


0.80 


0.74 


0.73 


0.14 


0.20 


0.35 


0.47 


impaired/enhanced 


eo- =0.27s-^ ^O.Ss^^ 


5.5 


4.9 


2.0 


0.35 


0.44 


0.45 


0.27 


0.23 


0.45 


0.48 


enhanced/impaired 


7ro- = 1.6s-^ ^2.5s-i 


5.4 


5.5 


0.35 


0.97 


0.67 


0.65 


0.22 


0.20 


0.40 


0.48 


impaired/enhanced 


WF-=0.65ifi ^l.Oifi 


5.3 


5.4 


0.85 


1.4 


0.59 


0.60 


0.21 


0.24 


0.42 


0.70 


-/enhanced 


VB- = 72 2S ^ 6.0 ifi 


5.8 


5.3 


2.1 


0.65 


0.63 


0.66 


0.43 


0.21 


0.52 


0.47 


enhanced/- 


iV_ = 6 -^ N- =5 


5.3 


4.4 


1.3 


0.37 


0.55 


0.58 


0.21 


0.22 


0.42 


0.46 


enhanced/impaired 



*The experimental values are taken from [12, 13, 14] as follows: The average stall forces for Wt II and Wt III are directly 
given in [12]; the other stall forces are read off from the diagrams in Fig. 3 in [13] and Fig. 2, 3 in [14] by applying the 
procedure described in the experimental procedures of [12]. The average run lengths are from Tab. 11 in [13] and Tab. I in 
[14], the average pause times from Tab. I in [13] and Tab. II in [14]. The average velocities for short and long runs have been 
read off from histograms in Fig. 7 in [13] and Fig. 4 in [14]. Missing values were not available. 
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Table 4. Mutation and regulation in lipid-droplet transport: results of the fit to the Drosophila lipid-droplet data that have 
not been used in the fitting procedure. In all phases and genetic backgrounds, the run length distributions can be fitted by a 
double exponential function, see Fig. 5A and Fig. 10 with the short and long decay lengths given here. Although not used 
for fitting, the simulation and experimental* values are of the same order of magnitude; the maximal error is 50%. 







decay lengths [/im 








short 


lor 


ig 






plus 


minus 


plus 


minus 


Wt II (sim.) 




0.099 


0.13 


0.95 


0.80 


Wt II (exp.) 




0.067 


0.098 


1.1 


1.1 


Dhc''-'"/+ (sim.) 




0.010 


0.13 


0.65 


0.61 


D/ic«-^7+ (exp.) 




0.088 


0.10 


0.78 


0.90 


Dhc^-^/+ (sim.) 




0.094 


0.083 


0.45 


0.33 


Dhc^'^/+ (exp.) 




0.074 


0.091 


0.40 


0.65 


Dhc^'-^/Dhc'-^" 


(sim.) 


0.059 


0.079 


0.31 


0.14 


Dhc^-^/Dhc''-^° 


(exp.) 


0.052 


0.044 


0.44 


0.21 


Wt III (sim.) 




0.13 


0.14 


0.48 


0.67 


Wt III (exp.) 




0.096 


0.083 


0.78 


1.1 



*The experimental decay length values are from Tab. I in [14], and Tab. II in [13]. 
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